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Abstract. Magnetic fields likely play a fundamental intermediary role between 
gravity and radiation in many astrophysical rotators. They can, among other 
things, 1) induce and be amplified by turbulence, 2) energize coronae, 3) launch 
and collimate outflows in "spring" or "fling" mechanisms The first is widely 
recognized to be important for angular momentum transport, but can also 
produce intrinsic variability and vorticity growth. The second leads to the 
production of high energy fiares, and also facilitates a test of general relativity 
from AGN observations. The third can operate from rotators with a large 
scale fields, though the origin of the requisite large scale fields is somewhat 
unresolved. I discuss these three points in more detail below, emphasizing 
some open questions. 



1. TURBULENCE, VARIABILITY & VORTICITY IN 

ACCRETION DISKS 

Accretion disks are a widely accepted paradigm to explain a variety of spectral 
features in sources such as active galactic nuclei (AGN), X-ray binaries, cataclysmic 
variables (CVs), dwarf novae and protostellar [1,2,3]. As gas orbits a central mas- 
sive source, internal dissipation drains the orbital energy, allowing material to move 
in and angular momentum out. Some fraction of the dissipated energy accounts for 
the observed luminosity. Micro-physical viscosities are often too inefficient so an 
enhanced transport mechanism, likely involving turbulence, essential. The obser- 
vational evidence for an enhanced/turbulent viscosity is least direct in the case of 
AGN disks, though it is natural to expect that the latter would be subject to the 
same turbulent viscosity generating mechanisms: As astrophysical disks are likely 
seeded with a magnetic field, the "Balbus-Hawley" or magneto-rotational instabil- 
ity (MRI) [3,4] ensues. This produces self-sustaining turbulence for flows with a 
radially decreasing angular speed that transports angular momentum in rotation- 
ally supported disks. Turbulence driven by the MRI transports angular momentum 



outward, unlike convection [3,5,6]. Rossby wave vortices may also transport angular 
momentum outward [7]. 

For low enough accretion rates, the dissipated energy may be primarily advected 
by hot protons in advection dominated accretion flows (ADAFs) rather than ra- 
diated, possibly explaining why some sources have engines that are surprisingly 
quiescent [9,10,11]. For optically thin, geometrically thick ADAFs, the disks may 
be physically similar to the coronae above the thin disks of Seyferts [12]. The thin 
disk+corona is something like a hamburger sandwich with the optically thin (geo- 
metrically thick) corona being the bun and the optically thick (geometrically thin) 
disk is the meat. The ADAF is a sandwich with no "meat." 

Alternatives to simple ADAFs suggest that quiescent luminosities might instead 
result from reducing the accretion rates into the inner regions [13,14,15]. Such 
models include advection dominated inflow outflow solutions (ADIOS) [13] and 
convection dominated accretion flows (CDAFs) ([14]; Narayan, these proc). These 
require a less than maximal viscosity [16]. 

Mean Field Disk Theory and Variability 

In general, the ubiquity of turbulence in disks has important implications for 
what the accretion disk equations mean. Analytic disk equations are mean held 
equations. While non-linear instabilities in thin disks have been extensively sim- 
ulated locally [17,18] a traditionally useful approach to disk models has been to 
swipe the details of the stress tensor into a turbulent viscosity of the form [19] 

vtb = aCsH ~ vthkb, (1) 

where H is the disk height, is the sound speed, Ith is the dominant turbulent scale, 
Vtb is the eddy speed at that scale, and ass < 1 is a constant. Use of this formalism 
requires a mean fleld theory, not a replacement of the molecular viscosity with the 
turbulent viscosity. Assumptions of azimuthal symmetry and steady radial inflow 
[1,11] require turbulent motions to be averaged over the time and/or spatial scales 
on which mean quantities vary. 

The required mean field approach is complementary to that employed in mean 
fleld magnetic dynamo theory, where the fleld is split into mean and fluctuating 
components, B = B + B', and the induction equation is averaged and solved. 
For the accretion disk case, the momentum equation must be similarly split, and 
the evolution equation for the mean velocity field derived. The result is that Uth 
represents a correlation of turbulent fiuctuations, that is Vth — (v'(i) • / v'{t')dt') 
[20,21,22]. 

Mean field theory has a limited precision, and thus predicts variability. Let the 
mean represent a global average over azimuth and half thickness, without averaging 
over the radius. Since all mean velocities are scalings of the Keplerian velocity, each 
radius is in principle "labelled" by its Keplerian speed. However, the averages over 
fiuctuations produce an uncertainty in this labelling. Since the luminosity at a 



given frequency depends on the radius of emission for accretion disk models, we 
can relate the uncertainty in the observed luminosity to the uncertainty in the 
radius. That is, 

AL^/L^ ~ \RdR{hnL^)\^R/ R = ^AR/R (2) 

where ^ can depend on temperature but ~ 1 for a range of frequencies for thin and 
thick disk models [23]. Let us estimate AR/R. The "error" associated with a single 

turbulent fluctuation of scale is reduced by A^^/^ for each averaged dimension, 
where N is the number of eddy spatial scales averaged over in that dimension. As 
observational data are taken in a time averaged sense, there is also a reduction that 
depends on observation time, tabs- The result is AR ~ ltb/[N^Nz{l + tobs/ttb)V^'^, 
where and A^^ are the number of dominant eddy scales in the and z directions, 
and ttb is the dominant correlation time. 

To make further progress, note that for fully developed MHD turbulence, near 
equipartition between kinetic and magnetic energy is generically reached, and is 
va ~ 2^/^i>t6 [3,17,18] in the case of accretion disks. The factor of \/2 comes 
from shear that adds a bit more amplification to the field. The MRI instability 
onset time scale is of order the dominant eddy turnover or cascade time and is 
Ub ~ kb/vtb ^ where fl is the rotation speed. Using these and (1) we have 

vy2n = assCsH. (3) 
Using Q,H ~ Cs for vertical pressure support, we then have 

ass ^ vl/2cl ~ ikb/Hf. (4) 

The next step is to note that A^^ = 27rR/{2lti,) = {TcR/HaH'^), where the extra 
2 on the bottom comes from eddy elongation in the direction from shear, and 
A^^ = H/2ltb = tt^/^^, where the 2 is from the 1/2 thickness. 
Collecting all of the above into (2) gives 

AL,/L, - ^ass{H/Rf'\l + ^Usy"^ ■ (5) 

For thick disk ADAF models, H/R ~ 1 and ass ~ 1 [8], large variability around the 
predicted luminosities can be expected unless tobs » When variability is not 
seen, or a systematic deviation from the theory is seen in a sample of observations 
at widely separated times or in different objects, the simplest ADAF type model 
may not be capturing the physics. Such systematic deviations seem to be evident 
in some of the large ellipticals [24] Other such quiescent accretor variations such as 
CDAFs (Narayan these proc.) or ADIOS type models [13] which involve outflows 
may be more appropriate there. 



Vorticity 



There are other effects of mean field theory in addition to the variabihty. If one 
considers the vertical disk averaging to be taken only over one hemisphere, then in 
addition to the scalar Shakura-Sunyaev turbulent viscosity transport term used in 
simple analj^ic accretion disk modeling, a pseudoscalar transport term also arises 
[22]. This term is analogous to that which appears in magnetic dynamo theory, 
and can lead to vorticity growth [21,25]. 

In the same way that the mean field magnetic dynamo characterizes an inverse 
cascade of magnetic helicity [26], the vorticity dynamo highlights some growth of 
vorticity on larger scales than the input turbulence. Enstrophy exhibits an inverse 
cascade in 2-D turbulence [27]. The growth of vorticity in primarily 2-D rotating 
fluids has been seen in nature (e.g. Jupiter [28]) as well as in simulation [29] and 
experiment [30]. Statistical mechanics approaches have modeled this [31]. Vorticity 
growth in sheared thin accretion disks has been studied less than th associated 
vortex evolution [32,33,34]. 

For an accretion disk in which mean quantities depend only on radius, Ref. [22] 
showed that 

dtUj = Vxckqw - V X {utbV X w) (6) 

in the frame comoving with the cartesian velocity tangent to the mean rotation, 
where the coefficients are 

ao = (r,/3)((a;'W-v'(°)) . . 

m = (rc/3) ( v'(o) • v'(o) + b'(o) • b'(o) ) ^ ' 

and b'^*^^ = B'^°)/47rp, with p as the mean density. (This equation presumes that 
1st order cross correlation terms vanish, that is (b'*^*^) ■ v'*^*^^) = (w''-'^-' ■ b'^^'') ~ 
• Vxb'("') = 0.) Note that Uth is not the only transport term. There is also 
the pseudoscalar ao term as in the magnetic field case. It is this pseudoscalar 
term which can lead to vorticity growth. Interestingly, the pseudoscalar for the 
vorticity, unlike for that of the mean magnetic field, has a kinetic helicity term 
without a current helicty term [22]. The can be parameterized as ao — qaggCs- 
One necessary condition for growth turns out to be g > H/ R. 

The simplest growth solutions [22] show a dominant growth scale ~ if, leading 
to intermediate scale vortices that should survive at least a vertical diffusion time, 
that is, ^ few x 1/a.ss orbits. For ags ~ 0.01, the resulting anti-cyclonic vortices 
may allow dust trapping, catalyzing planet formation when applied to star+planet 
system forming disks [31,33,35]. 

Note that this simplified model of intermediate scale vorticity growth cannot tell 
how many vortices grow, or where in height or azimuth these vortices are, only 
that there are growing solutions. This is because here the variables are averaged 
to depend only on radius. Note also that if the vertical averaging is taken over 



the full scale height, then the ckq coefficient should vanish because the psuedoscalar 
reverses sign across the mid-plane. Then vorticity growth could not be identified 
in this over simplified formalism. 

2. BUOYANCY, CORONAE & AGN IRON LINE 

PROFILES 

Coronae 

MHD turbulence likely involves spatial intermittency of the magnetic field [36], 
even in accretion disks. While the dynamics of intermittency is not fully under- 
stood, the random component of the field would be preferentially amplified at 
regions of strongest shear. Since small scale shear varies spatially and spectrally 
in a turbulent medium, the field strength would also be expected to to vary in 
correlation. The strongest magnetic field regions might form a kind of dynamic 
sponge, intermixed with "void" regions that have much weaker fields. The extreme 
situation, in which the magnetic field occupies a distinct volume from the thermal 
material, could maintain dynamical equilibrium when the ratio of the the average 
particle to magnetic pressure in the disk (= (3p) is large. To see this, consider the 
field to reside in magnetically dominated flux tubes. The pressure inside balances 
the external particle pressure. However, ii (3p » 1, the tubes would occupy only a 
small volume fiUing fraction (~ l//3p) [37]. Since the distance between tubes would 
be large, intersections with other tubes would be too infrequent to significantly 
load particles into the tubes. In contrast, if the magnetic volume filling fraction 
were large, frequent reconnection events over a large fraction of the tubes' longi- 
tudinal cross section would more easily mass load the tubes and lead to one phase 
medium. The amount of intermittency remains to be understood. Understanding 
this intermittency is important because extremely evacuated tubes rise at speeds of 
order their internal Alfven speeds and form coronae. Coronal magnetic dissipation 
is as important a problem for disks [38] as it is for the Sun [39] 

While coronae have been long thought to form above disks [38,40] analytic work 
often does not incorporate the turbulence and/or the origin of the initial large field. 
The usually invoked Parker instability favors wave modes which can be shredded 
by the Balbus-Hawley instability on time scales approximately equal to rise times. 
Simulations of turbulent disks show that coronae do form in turbulent disks [41] 
but the dominant mechanism, and its relationship to the process of formation in 
the solar corona is not fully sorted out. 

If the disk were not subject to MRIs and turbulence, the induction equation 
tells us that the disk field would grow linearly from shear. It would saturate at 
a value that is higher than the saturation value of the MRI. This is because the 
shear can amplify the field over a much longer time in a laminar disk; the field 
remains coherent longer. Field strengths of order the thermal pressure could incur 
before buoyancy (in this case by the Parker instability) drained the field into the 



corona. For MRI driven disks, the shear only operates on a single magnetic field 
filament for a correlation time (one rotation) after which that filament loses its 
identity. The field energy saturates by a factor of lower for the turbulent case. 
Also for the the laminar case, the corona would be more intermittent. This is 
because the buoyancy and the dissipation would likely get rid of the field faster 
than it builds up. For a laminar disk, the field growth would be linear in time 
whereas for the MRI disk, the growth is exponential. A laminar disk + corona 
should thus exhibit longer variability periods than a turbulent disk, but with larger 
amplitudes. The observational evidence in the case of AGN disks/corona seems to 
suggest a turbulent disk. While X-ray variabilities of factors of several are observed 
[42], factors of orders of magnitudes are not, suggesting a steady background level 
of coronal dissipation. 

AGN Iron Lines and Engine Geometries 

Scyfcrt galaxy X-ray spectra, are best modeled as a combination of direct emis- 
sion from a hot corona and reprocessed emission from a cold, optically thick ac- 
cretion disk [43]. The direct component results from inverse Compton scattering 
of thermal disk UV photons. The reprocessed component incurs as photons are 
scattered back onto the disk. 

If it weren't for the active coronae, we would lose a probe of strong gravity [44]. 
The reprocessed coronal emission includes the broad iron Ka fiuorescence fine of 
rest energy 6.4 keV, which carries information about the geometry and dynamics 
of the reprocessing material near the black hole [44]. ASCA has observed iron 
lines in ~ 18 Seyfert Is [45]. The best studied iron line is that of MCG-6-30-15 
[46,47]. In addition to the geometry, the iron line profiles are sensitive to the disk 
illumination law, the disk inclination, and the inner and outer radii [44]. Most 
work on reprocessing in AGN has invoked thin fiat disks with axially symmetric 
illumination laws representing a "point" X-ray source. 

Some AGN line profiles like MCG-6-30-15 are consistent with fiat disks [46,47] 
but the current data may be not precise enough to rule in or out non-axisymmetric 
engines in other cases. Such engines are also worth considering in order to provide 
robust comparisons to flat disk models. Ref [48] considered flnite disk thicknesses 
and Ref [49] considered concave disks. Warped disks have also been considered 
around Schwarzchild black holes [50]. There are plausible theoretical reasons for 
disk warping including radiation driven warping [51] and tidal warping [52]. On 
the observational side, water maser emission of NGC 4258 at O.lpc from the central 
engine (on larger scales than the inner accretion disk) traces a disk warp [53]. There 
may also be indirect evidence on these larger scales from observations of Seyfert Is 
which suggest that the broad line regions are not coplanar with the inner disk [54] . 
Dusty tori of Seyfert unification paradigms might also involve warped disks [55]. 

Warped disk studies reveal line features that are impossible for unobscured flat 
thin disks. First, shadowing of the source by the disk and shadowing of reprocessed 



emission by the disk blocks regions of the disk from contributing to the iron hne. 
Sharper red than blue cutoffs or very soft blue cutoffs can also arise. The latter 

characteristic is seen in some profiles of Ref [45]. The sharp red cutoffs result for 
large inclination angles, which is consistent with some Seyfert lis [56]. Second, 
non-axisymmetry of the disk means that line profiles can show time variability if 
the warp precesses around the disk. Third, there can be sharper peaks near the 
rest frequency compared to a fiat disk since concavity can offer more solid angle 
covering fraction. Fourth, apparent misalignment of central disk plane with the 
obscuring torus can be accounted for. 

For a concave disk, sharper peaks near the rest frequency can be accompanied 
by a total reprocessed emission fraction that is larger than 1/2 [49], where 1/2 is 
the maximum limit for a point source above a flat disk. This may play a role in 
ultra-soft narrow-line Seyferts (Brandt, private communication 1999), though there 
are other ways to achieve this enhanced reprocessing fraction. 

Line profiles from a distribution of dense clouds in an optically thin, geometrically 
thick disk may apply to ADAFs. Dense clouds formed from thermal instability can 
survive long enough to produce reprocessing signatures in otherwise optically thin 
flows [57]. This needs more investigation. 

3. SPRINGS, FLINGS & LARGE SCALE FIELDS 

Many large scale jets and winds in astrophysics including those of young stel- 
lar objects, microquasars, gamma-ray bursts, and AGN jets may be magnetically 
driven. Even supernovae may also involve magnetically driven bipolar outflows 
(Wheeler, these proc). How MHD jets work and where the requisite large scale 
magnetic flelds come from are integrated questions, but are usually studied inde- 
pendently. 

The large amount of work on MHD jet launching and coUimation will not be 
reviewed here (see [58,59]). However, note that the launching mechanisms could 
be divided into "spring" [60] and "fling" [61] mechanisms. In the former class, the 
jet is launched initially by toroidal magnetic field pressure. Imagine for example, 
a dipole magnetic field anchored in a star which incurs rapid differential rotation 
(such as the collapse from a white dwarf to a young neutron star as invoked in [62] 
for gamma-ray bursts) or in a supernova core collapse (Wheeler, Meier personal 
comm.). As the differential rotation winds up the field, the toroidal field pressure 
grows quadratically in time. When the pressure reaches some critical value, the 
field will act something like a coiled spring and can drive a strong torsional Alfven 
wave containing directed Poynting flux outward. This could in principle power 
a jet. Related mechanisms have been discussed for disks [60,63]. In this regard, 
note that in the case of AGN, it is not clear if the jet emission we see represents 
dissipation from instabilities at the edge of the jet, re-acceleration inside the jet 
along the bulk flow, or just emission from a very small number of particles carrying 
the currents [63] which support the magnetic fields. 



In the "fling" launch mechanism [61], the initial launch is driven by centrifugal 
force. The rigid field lines significantly weaken the effective gravitational potential 

when sufficiently inclined to the normal, allowing material fling out along poloidal 
fleld lines. Subsequently, before reaching the Alfven surface, the driving docs be- 
come magnetically driven as in the "spring" mechanism. Since simulation of such 
launching treats the base of the jet (e.g. its initial launch point in the corona) as a 
boundary condition, simulators often load the fleld lines with a little mass flow to 
get the process started [64]. The initial launching is different for spring and fling, 
but the ultimate collimation mechanisms could be the same, e.g. hoop stresses. 
The extent of the collimation appears to be sensitive to the boundary conditions 
of the outflow however [65] . 

MHD jet luminosity is fueled by the rotational energy. In systems which have 
both central compact rotators and disks, outflows could emanate from both e.g. 
[66]. In black hole systems, the relative contribution to the jet power from regions 
within and outside of the last stable orbit of the disk has been addressed [67] . Even 
if a Blandford-Znajek type mechanism is operating from the hole, a jet from the 
disk may in fact always dominate. However, the black hole spin can still influence 
the jet, since it determines the inner edge of the disk, and is ultimately important 
for understanding the magnctospheres [68]. 

Where do the required magnetic flelds [58] come from? The flrst possibility is 
that they are accreted. But this may not work for a turbulent accretion disk. 
Consider a turbulent disk threaded by a large scale vertical magnetic fleld. The 
fleld is subject to turbulent diffusion and may incur a net diffusion outward [69] 
(with some dependence on turbulent Prandtl number.) The role of reconnection 
may not yet be fully appreciated in this process. Without reconnection, the mean 
fleld is indeed subject to turbulent diffusion, but it cannot ultimately separate from 
the gas which has a systematic inward motion. In the absence of a topology change 
that can release fleld lines from the initial material they thread, the fleld would 
accrete on the diffusion time scale. 

If the flelds cannot be accreted then they would have to have been threading the 
central object before the disk formed, or be generated by a dynamo in the central 
object (if not a black hole) or disk. A traditional approach to the amplification 
of fields on scales larger than the scale of the input turbulence is the mean field 
dynamo [70] mentioned earlier. Only modes which have an initial seed field can 
be amplified. For an accretion disk dynamo, the limit of field energy density is 
the turbulent energy density. This limit is ass times the thermal energy density, 
which suggests that "spring" mechanisms are less likely from dynamo produced 
fields than "fling" mechanisms in turbulent disks. 

It is important to distinguish between "large" scale and "mean" flelds. Standard 
mean fleld theory is degenerate with respect to the topology of the fleld on scales 
smaller than the scale of the mean. Disconnected loops can have the same mean as 
a connected winding fleld line. In the formalism of mean field dynamo theory, re- 
connection is therefore not strictly required (though it is likely happening anyway). 
This is not commonly recognized. Said another way, neither the mean field nor the 



fluctuating component of the fleld are the topologically physical fleld. To generate 
jets, common wisdom holds that the mean fields actually do have to correspond to 
the topologically physical field. Perhaps this need not be the case if the Poynting 
fiux driving the jet is an average quantity: (E x B) = (e x b) + E x B, where 
the first term on the right, due to only fluctuating quantities, is usually ignored 
in this context (but see [71] where the coUimation is non- magnetic.) If we do ig- 
nore fluctuating components, then to magnetically launch and collimate jets by the 
"fling" mechanism, the mean flelds would need to be topologically large scale. One 
plausible way this could arise in a disk corona is if flux loops make their way to the 
surface and subsequent reconnection events inverse cascade smaller loops to larger 
loops [40]. 

The role of boundary conditions is particularly important for mean fleld dynamos. 
First, generating a net flux of the mean fleld in a quadrupole mode inside an object 
is accomplished by diffusion of the reverse flux through the boundary. Fast cycles 
(e.g solar cycle) of a dipole field also require boundary diffusion to change the flux 
inside. Incompressible simulations which employ periodic boundary conditions over 
the scale of the mean cannot see mean field growth because the induction equation 
then constrains the mean field to be time independent. 

In addition, ref. [72] showed that conservation of magnetic helicity means that 
the growth of large scale field with one sign of magnetic hehcity inside the system 
requires helicity of the opposite sign to diffuse out the boundary. Some results 
showing strong dynamo coefficient a quenching [73] may therefore highligh just an 
effect from the assumed boundary conditions rather than dynamical suppression. 
In general, to properly test large scale field formation in a turbulent disk one must 
really utilize a global study, with significant scale separation and diffusive boundary 
conditions. 



Closing Comment 

There is much to be learned on all of the above subjects by looking at the sun 
[74] , as others would also advocate [58] . Note that dynamos in the Sun may operate 
differently than in disks. For the Sun, strong shear amplification may take place 
below the actual turbulent zone, whereas in disks, the turbulent zone and the shear 
are the same region. However the sun has a large scale wind and with an active 
corona [74] , consistent with MHD outfiows along open field fines and x-ray activity 
resulting from dissipation of closed field lines. We should expect the same for a 
wide variety of turbulent astrophysical rotators. 
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